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Deconfinement transition in 2+1-dimensional SU(4) lattice gauge theory 
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A missing piece is added to the Svetitsky-Yaffe conjecture. The spin model in the same universality class as 
the (2 + l)d SU(4) theory, the 2d Ashkin- Teller model, has a line of continuously varying critical exponents. The 
exponents measured in the gauge theory correspond best to the Potts point on the Ashkin- Teller line. 
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1. FRAMEWORK 

The SU{4) Yang-Mills theory in (2 + 1) dimen- 
sions has a special status: its critical exponents 
at the finite temperature deconfinement transi- 
tion are not specified by the Svetitsky-Yaffe (S-Y) 
conjecture By measuring them, we can learn 
about the couplings of the equivalent spin model, 
i.e. of the effective Polyakov loop model. 

The S-Y conjecture says that, if the deconfine- 
ment transition of the gauge theory is second- 
order, then it is in the same universality class 
as the spin model having the symmetry-breaking 
pattern of the Polyakov loop. 

Assume that the (2-|-l)d S'[/(4) deconfinement 
transition is second-order - we will investigate 
this issue numerically next. Then the equivalent 
2d spin model has Z{A) spins, with orientation 
{0, tt} corresponding to the 4 Polyakov loop 
sectors and the 4 perturbative vacua of the gauge 
theory. A nearest-neighbour Hamiltonian, which 
gives the same critical properties as the 5^7(4) 
theory according to S-Y, will have 3 possible en- 
ergy levels for each link: Eq for parallel spins tt. 
El for perpendicular spins t^, and E2 for anti- 
parallel spins ti- Remarkably, the critical expo- 
nents of this spin system are known to vary con- 
tinuously with the ratio p = ^I-eI ^^"^ 
responding couplings of the Hamiltonian. There- 
fore, the S-Y conjecture does not tell us what the 
SU{4:) critical exponents are. Measuring them 
fixes the couplings of the effective Hamiltonian. 

A Z{4:) spin model with 3 energy levels per link 



is called a symmetric Ashkin- Teller model It 
is the symmetric case J — J' oi the Ashkin- Teller 
model, which describes 2 coupled Ising systems 
{(Ti, Ti} on a square lattice, with Hamiltonian: 
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(1) 



When J = J', it can be rewritten 



H = -2J ^ cos{e,-ej)~K cos{2{e,-9^)) 

<ij> <ij> 

(2) 

where now 9i — {0,±-|,7r} represents the orien- 
tation of an SU{4:) Polyakov loop. The ratio of 
excitation energies mentioned above, p = ^^ZeI ' 
is then equal to . Two special cases arise: 

• the 4-state Potts case: K = J ^ p = 0. No ad- 
ditional energy is needed to make an orthogonal 
spin pair anti-parallel. 

• the Ising^ case: K = ^ p ~ 1. It takes twice 
as much energy to flip one spin in an aligned pair, 
as to make it orthogonal. The name comes from 
the equivalent description as two decoupled Ising 
systems in Eq.([T|l. 

On physical grounds, one expects < p < 1, so 
that the two cases above are limiting cases. As p 
varies in the interval, the critical exponents also 
do, between the following limits: 
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where 7' and j3' are susceptibility and magnetiza- 
tion exponents for spins e^^^ . Since ^/v and (3/v 
are independent of p, these known exponents pro- 
vide crosschecks on our systematic errors. This is 
particularly valuable here, because of large, loga- 
rithmic finite-size corrections known to make the 
numerical extraction of Ashkin- Teller exponents 
quite challenging j^j- 

2. NUMERICAL STUDY 

We simulate the (2 -I- l)d SU{4) gauge theory 
using the Wilson plaquette action on a cubic grid 
of size 

X Nt, with Nt = 2, 3, 4 and L up to 40. 
To accelerate Monte Carlo evolution, we use as 
elementary update a mixture of pseudo-heatbath 
(in all 6 SU (2) subgroups) and overrelaxation in 
the full SU{4) group |4I5| . The latter requires a 
similar amount of work to 6 SU(2) overrelaxation 
steps, but gives a larger step size. As a result, 
the number of sweeps needed to decorrelate the 
Polyakov loop is reduced by a factor ~ 3. Up to 
10^ sweeps are performed on each volume. We 
analyze results for various couplings /3 near criti- 
cality together with multihistogram reweighting. 

This model gives a clear first-order transition 
on A't = 1 lattices. On Nt = 2, an old Monte 
Carlo study found a second-order deconfinement 
transition Unfortunately, the larger sizes 

which we simulated revealed that the transition 
is in fact first-order. Distortions of the plaque- 
tte distribution towards a double-peak structure 
appear in Fig. 1 for our largest sizes. Also, the 
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Figure 1. Distribution of the plaquette for various 
lattice sizes x 2. Deviations towards a double 
peak characteristic of a first-order transition start 
to appear at the largest sizes. 



r 




'/ 






y/v = 7/4 



Figure 2. Check of magnetic susceptibility ex- 
ponent ^ = |, without (left) and with (right) 
logarithmic corrections to scaling. 

Binder cumulant 1 — •^(^•lyi for the plaquette 
extrapolates as i — > oo to a value below 2/3. 
A likely explanation of these first-order transi- 
tions is the vicinity of a sharp, bulk crossover at 
/3 ^ 13.5, whereas the Nt = \ and 2 transitions 
occur at (i^ w 8.67 and 14.87 respectively. 

This forced us to consider Nt = i lattices, 
with correspondingly higher (i^- No double-peak 
structure is visible on the largest size considered 
(L = 32). However, the L — > oo extrapolation of 
the Binder cumulant still misses 2/3 by a small 
but somewhat significant amount. It may well 
be that the transition is still weakly first-order. 
Our current iVj = 4 results are consistent with a 
second-order transition, but do not reach as large 
volumes yet. Therefore, we present the critical 
exponents analyzed from the iVj = 3 data. 

The effect of logarithmic corrections to scaling 
can be seen in Fig. 2. They limit the usefulness 
of accurate, small-size data. 

The bulk crossover at /? 13.5 has a more per- 
nicious effect. It dominates the behaviour of the 
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Figure 3. — determined from the scaling of the 
specific heat, wihout (left) and with (right) sub- 
traction of smooth background caused by nearby 
crossover. 
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Figure 4. Correlation length exponent i/ deter- 
mined from the Polyakov loop Binder cumulant. 
The Potts value is strongly favored, even allowing 
for an error of ±0.01 in the determination of [3^- 
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Figure 5. Consistency of all determinations of Pc- 
V has the Potts value 2/3; corrections to scaling 
are fitted. 

specific heat over the singular piece, leading to 
an exponent a w (Fig. 3 left). We subtracted 
the bulk specific heat, measured on a 4^ lattice, 
to isolate the singular contribution. The expo- 
nent then becomes consistent with the Potts case 
(Fig. 3 right). 

All other exponents also favor the Potts case. 
Fig. 4 shows the scaling of where U is the 
Polyakov loop Binder cumulant. v remains con- 
sistent with the Potts value 2/3, even allowing 
for a large (> 2a') variation of (ic by ±0.01 and 
logarithmic finite-size corrections. 

Agreement among all observables is shown in 
Fig. 5. The pseudo-critical /3's obtained on vari- 
ous volumes all extrapolate to a common thermo- 
dynamic value /3c = 20.4356(41), with corrections 
of the form aL^~ (1+6/L), and v = 2/3, the Potts 
value. 



3. CONCLUSION 

The unexpected Nt — 2 first-order transition, 
together with expected but unwelcome logarith- 
mic finite-size corrections, have turned this simple 
problem into a numerically challenging one. 

Our results need to be made more precise, by 
simulating larger volumes. They also must be 
confirmed on a finer lattice, A^t > 4, where the 
transition presumably is second-order. In prac- 
tice, we can never exclude a weak first-order tran- 
sition. What we want is to reach physical volumes 
large enough to reliably determine effective criti- 
cal exponents, but small enough compared to the 
possibly finite correlation length at criticality. 

At this stage, the set of measured exponents 
favors the Potts case, which would be perhaps 
the most naive guess. 

Obviously, much remains to be done. The effort 
is worthwhile because of the insight gained from 
this simple case. Interest in Polyakov loop mod- 
els has been growing, because they can supple- 
ment dimensional reduction and provide an effec- 
tive description at temperatures nearer to Tc than 
the latter. The difficulty is to determine their 
many couplings, especially for large— A^ SU{N) 
theories. Lessons learnt here for SU (4) may help 
choose among the simplest SU{N) generaliza- 
tions, like A^— states Potts or Z{N) clock mod- 
els. In that respect, our preliminary result is 
perhaps surprising, because an effective A^-state 
Potts model for (3 + l)d SU{N) would lead to 
fc-string tensions independent of k at Tc. 
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